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Q/- Abstract 

The fundamental matrix solution of the quantum Knizhnik-Zamolodchikov equation 
associated with U q {sl2) is constructed for \q\ = 1. The formula for its determinant is 
given in terms of the double sine function. 



> 

on ■ 1 Introduction 
o 

The aim of this paper is to compute the determinant of the fundamental matrix solution of 
OO ■ — 
Q\ • the quantum Knizhnik-Zamolodchikov equation associated with U q (sl2) for \q\ — 1. 

In H and @, Tarasov and Varchenko studied the rational qKZ equation associated with 

the sh Yangian and the trigonometric qKZ equation associated with U q (sl2), respectively. 

The case studied in || is the limit q — > 1 of the case studied in J7|. However, there is a 

difference in these two cases. In the former, the unknown function ip(f3\, . . . ,/3 n ) is defined 

on C n , while in the latter the unknown function ip{z\, . . . , z n ) is defined on (C\{0}) n . If we 

set 

q = e p , z m = e p (1 < m < n), 

the limit q — > 1 corresponds to p —>■ oo. Therefore, the period pi in (3 m is lost in this limit. 

We will study the gKZ equation on C n with finite p. In [0, the solutions are single- 
valued on (C\{0}) n . The price for this restriction is that the multiplicative step p in z m of 
the difference equation is restricted by the condition < |p| < 1 (or \p\ > 1 in a different 
convention). We consider the case p = e p with real p and A. Namely, the additive step 
in j3 m is — Xi. Then, the solutions are not single- valued in z m . 

In the applications to physics, this difference corresponds to the difference of the models. 
In the application to the XXZ spin chain in the massive regime || , the deformation parameter 
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q satisfies — 1 < q < 0, and the multiplicative step is given by p = q~ A for the correlation 
functions and by p = q 4 for the form factors. On the other hand, in the application to the 
sine-Gordon model f|, the deformation parameter q satisfies \q\ = 1 and the additive step 
is given by A = 2n (i.e., p = q 4 ). 

We follow for the construction of solutions of the qKZ equations. We consider 

the evaluation modules of the U q (sl2) Verma modules with the spectral parameter z m and 
the highest weight 2A m . The unknown function ip(P±, . . . , /3 n ) takes a value in their tensor 
product. We will fix I and consider the subspace of weight 2(^^ =1 A m — I), which has 
dimensions d = f™ 4 ^ 1 ) . The qKZ equation is defined on this subspace. We will choose 
a basis of this space. The qKZ equation is then written as a system of equations for a 
d- dimensional unknown function. 

We will construct a non-degenerate d x d matrix such that each column solves the 
qKZ equation following a similar construction in ||. We call it the fundamental solution. 
The fundamental solution is constructed symmetrically with respect to the parameters p and 
A. Namely, each row solves another qKZ equation with 

— — 2ir/3m 

q = e A , z m = e A (1 < m < n), 

and —pi being the additive step in f3 m . We remark that in Q the co-existence of two systems 
of equations has already been discussed in the setting of elliptic equations. 

We give the fundamental solution in terms of /-dimensional integrals. To prove the non- 
degeneracy of the matrix we compute its determinant. The integrands of these integrals are 
meromorphic functions given explicitly in terms of the double sine function. We express the 
determinant itself as a simple product of double sine functions. In |||7j , similar determinant 
formulas are obtained. We follow their line of arguments in the evaluation of the determinant. 

We introduce another parameter p which couples to the s/2 generator h\ acting on the ra- 
th component of the tensor product. Therefore, the determinant is a function of the variables 
Pi, . . . , p n , Ai, ... , A n and p, A, p. The dependence on Pi, . . . ,/3 n is easily determined by 
using the two sets of qKZ equations. In order to determine the multiplicative factor which is 
independent of fti, . . . ,(3 n , we compute the asymptotics of the fundamental solution in the 
region P% « ■ ■ ■ « p n . This asymptotics is triangular, and the diagonal extries are the 
same kind of integrals with n — 1. 

The computation of the n = 1 integrals is done in two steps. The first step is to derive 
difference equations in p with steps — and The dependence on p is determined by these 
equations. Finally, the multiplicative factor which is independent of p is determined by 
taking the asymptotics when p — > zoo. 

The plan of the paper is as follows. Section 2 contains preliminaries of the qKZ equation. 
In Section 3, we construct the fundamental matrix solution. In Section 4, the asymptotics of 
the determinant is computed in terms of the n = 1 integrals. These integrals are computed 
in Section 5. The formula for the determinant is given in Section 6. A short Appendix is 
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given on the double sine function for the convenience of the reader. 



2 The qKZ equation 

Let q be a nonzero complex number which is not a root of unity. Consider the quantum 
affine algebra U q = U q (sl 2 ) with generators e^, q ±hi (i = 0, 1) and relations: 

q h iq hj = q h ]q h % ghiq-hi = q -hi q hi = 1? 
q hi ejq~ hi = q a *e h q K f 3 q^ = q~ aij fj, 

q ^i q 

fj] ®ij Zi J 

q — q 1 

e i e j ~ [ 3 ]q e l e j e i + [^]q e i e j e i ~ e j4 = ° (« ^ j), 

f!fj - *J?.I)L + ■ , > <J,f J ff - fjfi = (i ? j). 

Here a 00 = a n = 2,a 01 = a w = —2, and [n] q is the g-integer defined by [n] q = q q Z q -\ • We 
also use [n] q \ = Yll=ii m h- 
We use the coproduct 

A(q hi ) — q hi ® q hi , A(e,) = e t ® 1 + q hi <g> e i? A(/<) = /* ® g~ hl + 1 ® f. 

We denote by C/ 9 the subalgebra generated by ex, fx, q ±hl . Let A be a complex number 
and V A be the Verma module for U q with the highest weight q 2A . We denote the highest 
weight vector of V A by t/ ' and use the following basis {v^}k>o- 

exv {k) = [2A - k + l) q v^ l \ q hl v {k) = q ^~ k ) v ( k ) } f lV ( k ) = [k+ \} q v {k+l) . 

For a nonzero complex number z, we can give V A the U q module structure by setting 

eov = zfxv, q h °v = q~ hl v, f v = z~ l dv, (v G V A ). 

We denote by V A (z) the U q module defined in this way. 

Let Ai, A 2 be complex numbers, and v^^v^ the base vectors of V Ai , V^ 2 , respectively. 
For generic zx, z 2 , there exists a map R AlA Azx/ z 2 ) G End (V A (zx) <8> Vl 2 (z 2 )) such that 

R q AlM (zx/z 2 )vP ® 4 0) = ® 4 0) , (2.1) 
K lA2 {zx/z 2 ) o A(x) = a(A(x)) o R q AlA2 {zx/ z 2 ) for all a; G E^, (2.2) 

where <r is the permutation map a(a <g> 6) = 6 <g> a. Such a map R AlA _ 2 (zx/z 2 ) is uniquely 
determined by ( |2.1| ) and (2.2). 

Let V Ai ®V1 2 = ©z^o^a 3 l +a 2 -2« b e the decomposition of the U q module into the irreducible 
components. We denote by II; the projection to the component V A +A 2V Then we have 
(see, e.g., f§) 

K lA MM = HU«» X> II (2 - 3) 

«=0 i=0 y 1/2 



where 

RUM = £ S (1 n q2)2 i ^ 9 qhlh)K 

k=0 n,=i(i-^) 

Note that RjI iA A0) is triangular with respect to the basis (g> u^. 

Consider the tensor product W q = V^Zi)®- ■ -®Vl n (z n ). Let Rjife/zj) G EndW 5 be the 
operator acting on the i-th and j-th components as R\. A .(zi/ ' Zj) G End(V^.(zj) ® V\.( z j))- 
We denote by H m the operator on W q given by 

m— th 

H m = 1 ® • • • <g) hi <g> • • • ® 1. 
Let p, r be non-zero complex numbers. For m — 1, . . . , n, set 

^Oi,--- ,z n ;q,p,r) = R^ ^(pzm/ 2 m -i) • • • Rm,x(P z m/ zi)r 2Am ~ Hm 



X 



Rm,ni Z m/ Z n ) • ' • Rm,m+1 { z m/ z m+l) ■ 



The quantum Knizhnik-Zamolodchikov equation is the following system of equations for 
unknown function ip(zi, • • • , z„) that takes a value in W 9 : 

V^i, ■ • • ,pz m , ■ •• ,z n ) = K m (zi, ■ ■ ■ ,2»;g,p,r)^(«i, • ■ • ,2 n ), (1 < m < n) (2.4) 

In this paper, we consider the qKZ equation where the step p satisfies \p\ = 1. In this 
case, solutions to ( |2.4| ) are not necessarily single-valued with respect to the parameters 
Zi, . . . , z n , q,p, r. Therefore, we rewrite ( |2.4|) as follows. We set 



#l I - ) W=W ) ) - I 2f ) ) ) 

K m (fa,~-, Pn\ P. \ P) = K m {z[ P \ 4 P) ; q {p) ,P ip ' X) , r (A) ). (2.5) 



Then ( |2.4|) is equivalent to the following system of equations for ip: 

ip(Pi, ■ ■ ■ , P m - Xi, ■ ■ ■ , p n ) = K m (/3i, ■■■ , f3 n ; p, A, n)ip((3i, ■ ■ ■ , f3 n ), (1 < m < n). (2.6) 
For a non-negative integer /, set 

n 

Wf = {ve W q - h x v = 2(J2 A m - l)v}. 



m=l 



Since the matrix K m acts on for each I, the equation Q2.6|) for ip splits into the equations 
for the weight-/ component ipi £ ' . We denote by K m ^ the Wj 1 block of K m . Then, we 
have 

IpliPi, ■ • • , P m - M, • ■ • , fin) = K m ,l(Pl, ■ ■ , Pn\ P, X, (JL)l/)l{Pi, ••• , Pn), (1 < TTl < n).(2.7) 
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Let Wj 1 * be the dual space of with the coupling (v*,v) for v* G Wf* and v G . In 

n(p) ^(A) 

Section 3, we construct ^(p, . . . ,p n ; p, A, p) taking a value in Wy <8> Wy such that for 
any G Wf (A) *, 

^(A, . . . , A,) = • • • , p, A, /i)) g ^ (P) 

satisfies (|2.7|) and also, for any w* G wf P) , 

MPu--- , Pn) = (w*, %(P, . . . , (3 n] p, A, p)) G Wf ] 
satisfies (|2.7|) with p and A interchanged. 

3 The hypergeometric solutions 

For a non-negative integer Z, we set 

2p = i (hi ' ' ' Jn) ^ Z> ; Z m = / 

L m=l 

For L = (Zi,-- • ,/„) G 2^, set 

■^m = {h + ' ' ' + lm-l + 1, ■ • • , h + • • ■ + l m }, 

and define a map 7 L : {1, ■ ■ • , 1} — ► {1, ■ ■ • , n} by 

T, / -\ def . t-, r, 

We also define a partial order ^ on 2™ by 

L <Ll •<=>■ Z n < Z^, Z n _i + /„ < +l' n , . . . , Z2 + ■ ■ • + Z„ < Z2 + ' ' ' + C- 
We define a function Wj/ as follows: 
w[ p) (ai, • ■ ■ , an; fa, ■ ■ ■ ,p n ) = Skew (g { [\ai, ■ ■ ■ ,aj;^i, ■ • • ,A0) > ( 3 - 8 ) 
^W" ,a«;A,"- A) = ? wEm<m ' Ww II sh-K.-^-Tri) 



l<j<j'<l F 



x[J ei^-Vw+Vw") J] sh-^-^ + A^) J] sh -( aj - p m - A m7 ri) , . 

i=l V m<7 i (j) ^ m>7 L (i) ^ 

where Skew is the skew-symmetrization with respect to (ai, • • ■ , aO, he., 

Skew/(ai,--- ,a t ) = - ^(sgno-)/(a ffl , • • • , a CT J. 

' <res. 



We abbreviate w£ (ai, • • • , af, f3i, ■ ■ ■ , (3 n ) to Wj/((3\, ■ ■ ■ , /3 n ) (or w[ \a,\, ■ ■ ■ , «/)) when the 
dependence on the abbreviated variables is irrelevant. 
We set 



w L . 



7® =J2 Ct 

Lezp 

Now we define a pairing between and T^> . We use 

1 n 

<l(z; A ) = |<? ^, V _ A-n-V^r-^- A-rrV <P( X ) = \\ ~ fim] Kn) , 4> (x) = <f{x; -1) , 

(3.9) 



^2 («X — A7r)S'2(— — A7r) 



where ^(x) = ^(xjp, A) is the double sine function with periods p and A (see Appendix). 

In this paper, we assume that the auxiliary parameters, p, A, p and — A m (1 < m < n), 
are in R>o- 

Suppose that /, g are entire functions in the variables azi, . . . , a\. For fixed j3\, . . . , f3 n G R 
we set 

/(/, #) = f •■■ f dan-'-dan e ME ^ Qj JJ 0(aj) JJ - a f) 

Jc Jc j=i i<j<j'<i 

x /(«!, . . . , a I )flf(ati, ... , a { ). (3.10) 

Here the contour C for the variable atj is taken to be the real line R. Note that the poles of 
the integrand at 

(3 m - A m ni + piZ> + AiZ> , ay + m + p«Z> + AiZ> 

are above C and the poles at 

(3 m + A m iri - piZ> - AiZ> , acf -iri- piZ> - AiZ> 

are below C . As we will see in Section 4, there is a region of the parameters where the 
integral (ft.lO|) is absolutely convergent. 



The integral, in particular, defines the pairing I (w^ p \ w^) between G and 

Theorem 3.1. For G J 7 ^, we set 

^ = Hw[ p) ,w^)v { L p) where v[ p) = v[ h) <g> • • • ® u<W G V^Cz? ) ® ■ • • ® V^\z^). 
Lezp 

Then, ipi is a solution to fl2~~ 
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Proof. For L = (h, ... ,l n ), set L (m) = • • • , l m+1 , l m , ■ ■ ■ , l n ) and L = (l n , l u --- , Z n _ x ). It 
is easy to see that ipi is a solution to (|2.7|) if the following relations ( |3.11[) and (|3.12| ) hold. 



^ 41) (A? ' ' ' ) An+l) A? 



•Pn)v ( L P) 



(3.11) 



I{w { -f\f3 n - \i,/3i, ■ ■ ■ ,/?„_!) V A) (A, • • • ,Pn ~ Ai)) = e-^iofU 1 ^). (3-12) 



The relation ( |3.11[ ) follows from Lemma 5.2.2 in [[|. We prove ( |3.12j ) in the form 

I(g[ p \w^)(P u ...,p n + M) = e^ Xi I($\(3 n , ft, ■ • • , M K \Pi, • • • , AO). (3.13) 
First note that 



w 



{A) (ai, • • • , Oj + Ai, • • • , an) = (-l) n+l V A) («i, • • ■ , a h • ■ ■ , a{), 



W (A) (/?1,-- - tPn-uPn + Xl) = (-1)V A >(/V-- ,Pn-X,Pn)- 

From (16.59 )) we find 

<p(x — Xi; A) 



(3.14) 
(3.15) 



<p(x;A) 



sh -(x + Am) 
sh^(a; - A7rz - Ai) ' 



(3.16) 



Set 



i=i i<i<i'<* 



Then 



/tfU (A) )(A,-,W 



daa • ■■daif L (a 1 , ■■■ , a z ;/3i, • • • , /5 n )w (A) (ai, - - - , az;/3i, • • • ,/3 n ). 
First, we consider the case l n = 0. In this case, fi has no poles at 

OCj = P n + A n iri - piZ> . 

Therefore, when we make the analytic continuation (3 n — > P n + Xi, no poles of fz cross the 
contour C. From (|3.15|) and ( |3.16|) , we find 



h(Pl,--- ,P n + Xl)w^(P 1 ,--- ,P n + Xl) = f T (Pn,Pl,--- ,Pn-l)w^ (/?!,••• ,/3 n ). 



(A), 



Therefore, the equation ( |3.13| ) holds if l n = 0. 

In the case of l n > 0, if j G r^, the poles of fi at 

Oij = Pn + A n 7rz - piZ> 

may cross the contour C . In order to avoid this crossing, we shift the contour C to C + Xi 
for all j G r^. We note that /t has no poles at 

fflj = Pm~ Kniti + piZ>o (m = 1, • • • , n - 1), 
aj = aj/ + 7ri + piZ> (j' G T^, m = 1, • • • n — 1) 

for j G T^. Hence, this shift of the contours does not cause any crossing of poles. 

By changing the variables atj — > atj + Xi (j G T^) after this shift and using ( |3.14|) , ( |3.15| ), 
and ( |3.16| ), we get 



Jc Jc 

x h( a i~in+i, ■ ■ ■ Pn, pi, • • • ,/? n _i)w (A) (ai, • • • ,aj;/3i,-.. ,/3 n )(3.17) 



Symmetrizing the integrand of (|3.17| ), we obtain ( |3.13| ). □ 



We define the fundamental matrix solution by 

. . . , p, X, A = J (4 P) ; w$)v%> ® «W G ® wf >. (3.18) 

L,L'£Z™ 

Then, it has the property announced in Section 2. 

4 Asymptotics of the Determinant in fii, • • • , /?„ 

In the following sections we calculate the determinant 

L,L'eZV 



JL sh^(/3 m -/5 fc + (A m + A fc -j)«) N 
X sh z(p m -f3 k - (A m + A fc - j)7ri) ^ 



-1 



fe=m+l 
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(4.19) 



A(/3i,--- ,A0 = det [/(w^,^) 

From Theorem 3.1, we find 

Di(Pi,~' ,Pm-Xi, ■■■ ,p n ) 

— — — — — = detK mi ,(Pi, • • • , Pn! P, A, p) 

U l\Pli " ' " j Pmj " - - , Pn) 

A( n^" ,/3 V PV 'Vf ra) = ^ K ™APn ■■■ ,Pn',X, p, p). (4.20) 
Using the formula fl2.3|) , we see that 

detK mi i(fti, ■■ ■ ,/3 n ;p, X,p) 

[ ' f} V fc =i sh ^ - p k - Xi - (A m + A, - j)m) 



where (?) is the usual binomial coefficient. 
We consider the following function 



11 + 

ii LJ<L_ n s *w™ - m - (a™ + Aw - * • J 



By using ( |6.59| ), we can check that Ei(/3\, - ■ ■ , f3 n ) satisfies both (|4.19|) and ( |4.20| ). Therefore, 
we have 

Proposition 4.1. 

•••,/?„)= Q(p, A, //; A 1; • • • , K)Ei{Px, •••,/?»), (4-22) 
where cj(p, A, p; A 1; • • • , A n ) is a constant independent of Pi, . . . , ft n . 

In order to determine ci(p, X, \i; Ai, • • • , A n ), we consider the asymptotics of Di/Ei as 

A,-AeR, A < (4.23) 
Hereafter we use the notation ~ as follows: 

f(Pi,--- ,Pn) ~9(Pi,~'Pn) & f( %''"'%i ^ 1 in the limit. 
We use the abbreviation P mm > = P m — P m ' ■ From (|6.57| ) , we find 

Eitfi, •••,&)- n p% jt (4.24) 



where 



* — (c + r ot-^s*- - (3 ("»^ o))- <4 - 25) 

27T 2 ( n \ n 

Pl{P\i ■ ■ ■ , Pn) = exp (~ T~ [ ^mA m '|/3 mm '| — Imlm'Pm'm J ~ H ^mA) 



\m,m'=l m<m' / m=l 



Here we used the equalities 

E( n + I — 1\ . 
l mi ---lm k =[ +k _i) V-<mi<-- - <m k <n). 



We have 

\p L l[w { [} ,wf) . (4.26) 



A(A,---,/9n) l det 



Ei(Pi,--- ,P n ) ~ dj 



In the following, we consider the asymptotics of PlI(wjj ,wj* ). We note that 
I(w$\wp) = ^(sgn a)I(g$(a ai ,... , a ai ), #[ A) (ai, . . . ,«/)). 



By changing the integral variables a p — > a p + /3 7 £( p ), we have 

1(9$ (<x n > ■■■ , a<n), #i A) (ai, • • • , aj)) = / • / d«i • • • da^U"^^ 

Jc Jc 

i 

j=l l<j<j'<l 

x ^(«ai + /V(<n)> • • • + ^ L {n))9L i a i + Ay £ (i)> • • • + /V(o)- 
We set 

J£ )L ,(ai, . . . ,a t ;p 1 ,... , (3 n ) = P L (/?i, • • • ,/?„) x (the integrand of (gg^)). 
From QS.58| ), we have 

(s — ► oo, ±Rex > 0) 



(4.27) 



p + A + 2Att 
AJ ~ exp I =F7r ; a; 



pX 



Therefore, if llmxl < K, we have an estimate 



\(p(x; A) | < 7^exp ( -tt 



p + A + 2A7r , 



\x 



(4.28) 



(4.29) 



where 7# is a constant independent of x. Set £(2) = a; + \x\. Using (|3.9|) , and Q4.29Q , we 
obtain the following uniform estimate in the asymptotic region (|4.23 ) . 



|J£ L ,(ai, . . . ,oti; ,/3 n )| < 7exp(pJ^ ^ + — ^ |a 3 - - 



2tt 5 



/ 



a f \ 



3<f 



(4.30) 



Here 7 is a constant independent of a%, . . . , aci, fli, . . . , (3 n . 

Now, it is easy to see that J£ L , is uniformly integrable in ( [4.23D if 



2tt2 



2tt 2tt 2 



m=l 



(4.31) 



m=l 



For simplicity, in the following calculation, we assume that p and A are sufficiently large. 
We have, in particular, that the region of convergence (|4.31 ) is not void. 

If L -fs L', then for any a G Si there exists j such that 7 i (<7j) > 7 L '(j)- Therefore, by 
Lebesgue's convergence theorem, we have 
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in the limit ( |4.23| ). Therefore, we get 

A(#L> • • • ,Pn) 1 TT D T ( (p) (A)x 



(4.32) 



(p) „„( A h 



Let us calculate the asymptotics of PlI{w£ ,w l ). Note that 



7> P ) < y (P) f o r all p <{=► a E S L = S h x-xS^C^. 
Therefore, we have 

p l i{w { 1\w { l ] ) ~ ^ E ( s s n(T ) y* ' *y dai ■■■ da i ^l,l("i> • • • A, • • • > AO- 

Using fl3~9l) , and ( |4.28|) , we can calculate the limit of J[ L (<r 6 S (l x ••• x S in ) as (|4.23|) : 



{p) (A) / 1 !.../J(g(P)gW)E m<m ^ m ^ 



4 i («- 1 ) + El< m <m'<„ *mV /! 



X 



2tt 



't71— 1 



m=l 



j=m+l 



P + X 

PA 



vr 1.(4.33) 



Here F£(x) is given by 



F A (a;) 



/ • - - / doti • • • da k e xT - k ^ Qj JJ A) JJ ^(ot 3 - - ay) 

i =1 i<j<j'<k 

x Skew I jQ sh — (oty — otj — iri) J f | sh T( a j' — a i ~~ 7ri)(4.34) 

\l<j<j'<k P ) Kj<j'<k 



From flOg ), fl43g ), and flOSD , we get 



Proposition 4.2. 



( g (p) ? (A))(S)( n ^ 1 )+("^ 1 ) EIU A m i 

c z (p,A,/i;Ai,--- ,A n ) = — 1 ^„ + ,_ ni ^ Vw+ ,_ n „ - ^„ + ,_ 1 v ||(j!) ? 



(n + l-j-2\ 



2tt 



n n^^ia 



4 »(..-i)(»r l )+fi)(",t;: 1 )(fl)(",S') y 

2 /m— 1 n 

E & - A ^ 



LG.Z," 771=1 



vi=i 



J =771+1 



p + A 
pA ' 



-7T 



(4.35) 



where Fj^ 1 (x) is defined by ( |4.34[ ) 



5 Calculation of the integral 

In this section, we find an explicit formula of F^x) 
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The integral ( |4.34| ) is absolutely convergent if 

iRexl < - + ^-^-(/-l-A). (5.36) 
p A pA 

In fact, F t A (x) is analytically continued in x to the whole complex plane. To see this, we 
derive difference equations satisfied by F A (x). Namely, we relate the values of F A at x + — 
andi + j to i^ A (x). 

Proposition 5.1. The function F^x) satisfies 

F t^-D l [ ch(§x + ^(k-A))' (5 ' 37) 

g^+j) _^ ch(fx-f(fe-A)) 

^-f) Mch(fx + ^-A))- ( ' } 

Proof. Because of the symmetry between p and A, it is enough to show ( |5.37 ). 
We set 

i 

f A (x\a u --- ,a t ) =e x ^U^ ^ <p(a i; A) JJ ^(ay - a?) 

i=i i<j<j'<i 

x Skew! ] [ sh — (atf — ctj — iri) J, (5.39) 
Vi<i<i'<; ^ / 



7T , 

x yCtf 

l<j<j'<l 



h(aci, •••,«/)= J J sh -r{oi,> — n , ■ — tt/ ) 



The poles of / A in the variable ctj that are lying above the contour C are at — Ani + piZ>o + 
AiZ> and ay + iri + piZ>o + AiZ> . 
Consider 



-ffc = /''■/( / ~ ) ' ' ' <^ a «/ A ( a; | a l; 

ic ic \Jc+pi JcJ 



x h(a u - • ■ ,ai)sh?-(ai + A-Ki)e * e f'" ; ^=1°^. (5.40) 
A 

The integrand has no poles inside the strip < Ima/ < p, and the integral is absolutely 
convergent if |Resc| is sufficiently small. Therefore, we have 

4 = 0. 

Now we transform the variable ai to ai + pi so that the contour C + pi is modified to C . 
From ( p.9| ) and ( |3.16| ), we have 

/ A (s|a»i,-- - ,q f + /rc) = ^ shf(a;-A7rc) ^ sh | (a,- - a, - ttz) ^ 
/ A (x|ai, • • - , ct{) ' sh |(o;/ + A7rz + pz) ^ sh \{cti — ctj — ni + pi) 
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Therefore, the equality I k — gives rise to 

0= /■■■/ da l ---daif k {x\a l ,--- ,ai)e--^^^-^ a ^ 
Jc Jc 

^_iy e *pi s h - Km)h{a h an, ■ ■ ■ , aj_i) - sh -(a* + K-Ki)h{a u ■■■ , aj) J (5.42) 



x 

We use the notation 



2-7T A7T ? /\\ 

a j — e x > r = e * , g = g l j . 



Symmetrizing the integrand of ( |5.42j ), we get 



(5.43) 



= y • • J dai - ■■daif A (x - j\<*i, • ■ ■ ,ai)A k (x\a u ■ ■ ■ ,a t ), 
where 

A k (x\ai, ■■■ ,ai) 

= Skew {(-e x ' p? (r _1 ai - r)a k+1 • • - a; - (raj - r _1 )a fc • • -aj_i) /i(a 1( • • • ,«/)}. 
Since (a^ — g _2 a J+ i)/i(o!i, • • • , a{) is symmetric with respect to aj and Oj+i, we have 

Skew (djh(ai, ■ ■ ■ , a/)) = g _2 Skew (a 3 -+i/i(ai, ■ ■ ■ , a;)) . (5.44) 
Using ( |5.44|) repeatedly, we get 



Skew (aia k+l a k+2 ■ ■ ■ aih{a u ■■■ , a*)) = g 2(fc 1} Skew (a k a k+ i ■ ■ ■ aih(a t , ■■■ , an)) , 
Skew (a k ■ ■ ■ a^h^, ■ • • , on)) = g~ 2(i ~ fc) Skew (a k+1 ■ ■ ■ aih(a u ■■■ , an)) , 

and therefore 

A k (x\ai, •••,«/) = -(r^V'^g -2 ^ -1 ) + r)Skew (a fe a fc+ i • • • a t h(ai, ■■ ■ , a*)) 

+ (re** + r- 1 g- 2 ^- fe ))Skew (a fe+1 • • • a//i(a x , • • • , a,)) . (5.45) 

Substituting ( |5.45| ) for ( |5.43|) , we get 



t l e xpi q 2(fc 1} + r) / ••• / da x ■ ■ ■ daJ A (x - -|a)Skew (a k ■ ■ ■ a t h(ai, ■ ■ ■ , on)) 
Jc Jc * 



{re xpi + T-\- 2{l ~ k) ) 



J • ■ J dai • ■ ■ daif A (x - ^|«)Skew (a k+1 ■ ■ ■ aih(a x , ■ ■ ■ ,aj)) , 

(5.46) 



and therefore, 

i-i 



Y](r- 1 e xpl q- 2k + r) / ■•• /dai • ■ • donf A {x - ^|a)Skew (a x • • • a t h{a 1} ■ ■ ■ ,a,)) 
k=o J c Jc * 

= Y[(re xpi + r- 1 q- 2k ) I ••■ f d ai ■ ■ ■ danf A (x - j\a)Skew (h( ai , ■ ■ ■ ,a,)). 
fc=0 c ^ 

If | Rex | is small, both x ± ? lie in the region ( 5.36| ). Therefore, we obtain ( p. 371 ). D 
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Consider 



Gf(x) = TT — — t 5— — t r-, (5.47) 



where ^(x) = ^(xjp, A). We can check that Gf also satisfies ( |5.37| ) and ( |5.38| ). Therefore, 
we obtain 

Proposition 5.2. 

F l A (x)=c l (p,\;A)G?(x), (5.48) 
where q(p, A; A) «s independent of x. 

In order to determine q(p, A; A), we consider the asymptotics of F A and Gf as x — > +200. 
From (|6.58 ), we have 



k=0 J 



G\ (x) ~ exp [ — xTii (A — fe) ) , (x — > +zoo). 

Therefore, 

q(p, A; A) = lim exp ( xm V(A - k) ) F^x). (5.49) 



1-1 



fc=0 

We calculate the limit ( |5.49|) . We use the following equalities. 



Skew ( Yl sh^(a f - aj -m) \ = -^p JJ sh 

\l<3<j'<i P / ' Kj<j'<l 



TT , 

a.ji — aj), (5.50) 



P 

\1<3<3'<1 / i<i<i'<i 

V 3=2 P 3=lj'=j+2 P 

- Y[ sh-(a f -aj). (5.51) 
' i<i<j'<« P 



Combining Q4.34j ), ( |5.50|) , and ( |5.51| ), we have 



z-i 



exp xvri^(A - k) \F( 



k=0 



q^l- i-»(i-l)(J_2A-2) 

/! 



x = — — e 2 p 



/ • • • / dai • • • da t e x ^U °y+*<£&(A-k)) JJ . ; A ) JJ ^ _ a/ ) 

l 1-2 i 

x e~ p ^ J=1 aj Y\ sh — («j + Am) JJ j | sh — (a,-/ — aj — iri) ] J sh — (aj> — aj). 

3=2 P j=lj'=j+2 P l<j<j'<l 



(5.52) 
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We denote by I- the integrand of ( 5.52[) . The integral is zero in the limit x — > +ioo if 



/ i 1-1 \ 

Im a i + ni 5Z( A ~~ k n >0 - 

\j=l k=0 J 



(5.53) 



Therefore, we can compute the integral by shifting the contour in the positive imaginary 
direction, and taking the residues at the poles which we are crossing. Note that we assume 
that p and A are large. Therefore, the only pole we cross is the one at a\ = Km. Thus we 
get 



lim 

x— *ioo 



C JC 



da\ ■ ■ ■ dairf = 2m' lim 



C-Ani+iO JC-Am+iO 



da 2 ■ ■ ■ daiies ai= _AnJt- 



(5.54) 



Now the integrand has a pole at a% — (1 — K)m. Repeating a similar argument, we obtain 
lim / •■■/ dati ■ ■ -doiitf = (2iri) l res ai= n-i-A)m • • • res ai= _A7r»^ A - (5.55) 

From ( gg|) , (|532|) , ( |5~55D flOftJD and (|6T6T| ), we get 
Proposition 5.3. 



ci(p, A; A) 
where S^x) = S^xjp, A). 



(a); 



n 



S , 2 (7t)a/pA 



4(2)/! f- = 1 1 ^ 2 (A;7r)5 2 ((A;-2A-l)7r) ! 



6 Formula for the determinant 



Combining (fO|), ( ggg ), flOSD , QQ7D , (|5^8|) , and (|5T56|) , we obtain the following formula 
for the determinant A(/?i, • • • , AO- 



Theorem 6.1. 



det 



L,L'&Z 



( n\ ( n-f-Z — 1\ 1 / n+i — 1\ n * 
( ? (p) g (A))U)(n +1 ) + ( „ i^-^^^v^A) 



X 



I n 

nn 



4«("- 1 )(" + n" 1 ) + ( n r)C^ 1 )(/!)("^ 1 ) 



n + i-j-2 N i 
n-2 J 



i™=i Vnti S 2 (kK)S 2 ((k - 2A m - 1)tt) 



2-1 

n 

j=0 



^(^-(E^ = lAm-J» 

^(^ + (E^=iA m -j> 



("n-T 1 ) 



A)- 
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Proof. It remains only to calculate q(p, A, /x; A 1; • • ■ , A m ) by using ( |4.35| ), ( |4.21| ) and ( |5.56j ). 

It is easy to see that 

n 

n n^ A i A -) 



LeZJ 1 m=l 



n+l-l\ , . „ (n+l-j-2\ 

In/ r .-i , r .-i , \ V n-2 J 



4K »+i J jla ^-j n , =i 5 2 (/tt)S 2 ((/ - 2A m - 1)tt) 

From (|5^7| ), ( |639| ), and ( ggg ), we have 



G t: (a* + ^ Eft - - E & - j - ^ 



5 2 (g/i-7r(A( m )-/ m -/)) 



TT 



where A( m ) = £™ =1 A, - £™ =m+1 A J; and Z m = £JL i Z 4 - EJ= m +i 
Now we rewrite 

\ X\ \\ Y — Z — . 6.56 

L% i=i to - (A(-D - l m + k)rr) 

Let us consider the following sets: 

T± m) = U {T m ±k\k = 0,1,- ■■ ,l m -l] (m=l,...,n), 

Lez» 

where |J means a disjoint union. For a G Z, we set 

multi m) (a) = #{teTi m) |t = a}. 
We set mult± (j) = mult^ i+1 ' ) (j) = 0. Then, we have 

p\ xmultV m) 0')-multL m+1) (i) 



ra= nn^(^-(A (m) -j> 



m=0 jeZ 



M V^(^-(A (0) +J>) 



("n-l 1 ) 



by using 



mult+^j) = mult^(— j) 



(^T 1 ), if J=0,1,---,Z-1; 
0, otherwise, 



mult^O") - mult^ +1) (j) =0, (m = 1, ■ ■ ■ , n - 1), for all j E Z. 
This completes the proof. □ 
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Appendix Double sine function 

Here we summarize the property of the double sine function S 2 (x) = S 2 (x|(^i , 002) following 

@- 

We assume Reu^ > 0,Reu 2 > 0. ^(x^i, u; 2 ) is a meromorpic function of x and sym- 
metric with respect to U\,uj<2. Its zeros and poles are given by 

zeros at x = co>iZ< + u 2 Z< , poles at x = u>iZ>i + 0J?7'>i- 

Its asymptotic behavior is as follows: 

logS-^x) = 

= ±m ( — - ^±^x --f^ + ^ + 3^) (x-oo,±Imz> 0).(6.57) 

\2cuiuj 2 2uilu 2 12 \a>2 ct>i / / 

This implies 

log £2(0 + x)S 2 (a — x) = ±ni -x + o(l), (x — > 00, ± Imx > 0). (6.58) 

UJ1U2 



The following formulae hold: 

S 2 (x + UJi) 



(6.59) 



S^x) 2 sin 22' 
S 2 (x)S 2 (—x) = — 4sin — sin— , (6.60) 

S 2 ( x ) = — ^^x + Q(x 2 ) (x->0). (6.61) 
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